Transport signatures of Floquet Majorana fermions in driven topological 

superconductors 



Arijit Kundu and Babak Seradjeh 
Department of Physics, Indiana University, 727 East Third Street, Bloomington, IN 47405-7105 USA 

Floquet Majorana fermions are steady states of equal superposition of electrons and holes in 
a periodically driven superconductor. We study the experimental signatures of Floquet Majorana 
fermions in transport measurements and show, both analytically and numerically, that their presence 
is signaled by a quantized conductance sum rule over discrete values of lead bias differing by multiple 
absorption or emission energies at drive frequency. We also study the effects of static disorder and 
find that the quantized sum rule is robust against weak disorder. Thus, we offer a unique way to 
identify the topological signatures of Floquet Majorana fermions. 



Introduction. — The nonlocal quantum order character- 
izing the topological state of a gapped medium often ne- 
cessitates the existence of topologically protected gapless 
states bound to bulk defects or the edge with a topologi- 
cally trivial medium where the gap closes. The detection 
of these topological bound states is, therefore, a primary 
probe of the topological state. It has recently been under- 
stood that topological bound states may arise as steady 
states when a topologically trivial system is driven pe- 
riodically [2 m . In the superconducting state, these are 
equal superpositions of electrons and holes known as Flo- 
quet Majorana fermions [3l [4] that exhibit non-Abelian 
statistics [3 [6] and can be used for topological quantum 
computation [7]. This possibility expands the systems 
and conditions that realize Majorana fermions as emer- 
gent quasiparticles [5HTB]. but also poses fundamental 
questions as to how to detect and possibly manipulate 
such steady states. In particular, since the driven sys- 
tem is not in equilibrium, the experiments probing the 
equilibrium response of static Majorana fermions |17H20| 
cannot be used directly for this purpose. 

In this paper, we address these questions by studying 
the non-equilibrium transport properties of Floquet Ma- 
jorana fermions. We show, both analytically and numer- 
ically, that there is a quantized conductance sum rule, 
which we dub the "Floquet sum rule," whenever Floquet 
Majorana fermions exist. The Floquet sum rule naturally 
generalizes the quantized zero-bias conductance of static 
Majorana fermions [2TH26] . Moreover, we show that the 
Floquet sum rule is robust against moderate static disor- 
der, owing to its topological character, while other peaks 
get suppressed. Remarkably, this suggests that disor- 
der, usually detrimental to electronic properties, can be 
used as a "sieve" to find Floquet Majorana fermions. 
Transport studies in irradiated graphene, where the Flo- 
quet topological insulator was first proposed to exist [1], 
suggested quantized transport in the driven system is 
possible in certain geometries and for large drive fre- 
quencies jMJ [37] . We use a systematic Green's function 
method that extends the previous studies to supercon- 
ducting systems in any frequency range, and can, in prin- 
ciple, incorporate the effects of interactions. 



Though our results are applicable to any realization of 
Floquet Majorana fermions, systems of cold atoms could 
prove specially useful in this regard due to a high de- 
gree of design control and newly developed experimental 
probes of their dynamics, such as single-atom imaging, 
tunneling, and transport |27H33j . Disorder can be in- 
troduced in cold atom systems controUably |34j |35] and 
could, therefore, play a key role in the detection and ma- 
nipulation of Floquet Majorana fermions. 

Model. — We study the model Hamiltonian H{t) — 
Hw{t) + Hc + H\, where the last term describes the leads, 

HUt) = ll'A{th, (1) 

is the Hamiltonian of the system (wire) in the Ma- 
jorana basis 7''' — (71,- •• ,727v) with a real, skew- 
symmetric matrix A{t), and the contact Hamiltonian 
He = + h.c. with a^^ is the row of electronic 

creation operators in lead A, and a contact matrix. 

Our analytical results are presented for a general re- 
alization of Majorana fermions, but we do have to pick 
a specific model for numerical calculations. Perhaps the 
simplest model system that realizes emergent Majorana 
fermions is a quantum wire with superconducting pair- 
ing in a spin-polarized electronic band O [38] . This sim- 
ple model can be realized in solid state [T31 [T71 [TB] 
and potentially in cold atom systems [31 [35l442|. Mul- 
tiple bands in the wire can be included in 7 such that 
N = Ni,L, where Nt, is the number of bands and L is 
the number of sites. Each band furnishes two Majorana 
operators (7ri7 7r2) = site r. The contact matrix 

is K"^ oc (a'^c^) (E) (1, i) in the Majorana basis, where the 
second factor acts on the Majorana indices 1 and 2 at 
each site. 

For the single band wire {Ni, — I) the only nonzero 
elements of A are 

Ar.r — — iP'r^yi ^r,r+l — ^r^x -\- iw^-Oy, (2) 

where the real parameters ^r, and Wr are, respec- 
tively, the chemical potential at site r, the supercon- 
ducting pairing, and the hopping integral on (r, r + 1) 
bond, and (ct^;, Cj,, CTz) are Pauli matrices. The static. 
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L = 150,A/i« = 0.6,/Ji/U/ = 1.25,/i2/W = 0.55,(i/f2 = 1.0 




FIG. 1. The evolution of the Floquet spectrum of a quan- 
tum wire with L sites for a square-wave chemical potential, 
with (a) Q,/W and (b) ii\/W. The arrows in (a) show the 
frequencies used for the transport calculations in Fig. |3] 

uniform wire with bandwidth W — 4|Rew| has a topo- 
logical phase transition at 2|/i| = W where the gap closes. 
There is an unpaired Majorana fermion at each end of 
the wire with zero energy (the energy is referenced to 
chemical potential of the wire) in the topological phase 
2\ii\ < W and none in the trivial phase 2\^\ > [Sj. 

When the system is driven with period T = 27r/f2, 
the general solution of the time-periodic Schrodinger 
equation (we set h = 1) H{t)\^{t)) = idt \i^{t)) is 
found in terms of the Floquet functions Itpait)) = 
e~"°* |0Q(t)), where \cj)a{t + T)) — \<l)a{t)} is an eigen- 
ket of the effective Hamiltonian H^sit) = H(t) — idt, 
Hcs{t) \4'a{t)) = Ca \4>a{t))- The quasicnergics are 
restricted to (—17/2,17/2] by the map + kQ, 

|(^„(t)) e'*^"* |'/>a(i))- We shall compute the Floquet 
spectrum {e} using the evolution operator U{t) |V'(0)) ~ 
\ip{t)) and constructing the Floquet Hamiltonian Hp = 
(VT)log[C/(T)]. Then, Hf\MO)) = |0a(O)). The 
periodic eigenkets can be resolved in a Fourier series 
= Efc e-^'^"* l^e^)). We shall use a shorthand |0„)> 
for vectors in the extended Hilbert space spanned by 
with the inner product = J^k {(t>'^^^\<l)^^'') = 

(^'(010(0) dt/r[B]. 

Floquet Majorana fermions. — Floquet Majorana 
fermions are bound states with quasienergy cq — 
or = fl/2 3J. The particle-hole symmetry, 

ffj, = — i?w; requires Hp — —Hp, so the quasiener- 
gies come in pairs (eq,— ea). In the Nambu basis. 



-73 ^ 1 \^ ^ 1;" ) ' ^^"^ ^^^"^ Floquet 
Majorana fermions vo{t) = u5(<) and v.^{t) — e'^*u* (t). 

In Fig. [l] we show the Floquet spectrum of a wire 
with a square- wave periodic chemical potential, fJ-{t), al- 
ternating with frequency between /ii and /i2, respec- 
tively, over time intervals ti and ^2 in each period. Note 
that ^{t) is not in the topological range at any time. As 
il/W decreases, the Floquet band spreads and e crosses 
Ett, giving rise to a pair of Floquet Majorana fermion 
bound states. Additional gap-closing level crossings lead 
to the appearance or disappearance of sq and e^^ Floquet 
Majorana fermions. 

One expects that Floquet Majorana fermions appear 
for 17 < where exchange processes with energy 17 
become important and result in qualitative differences 
between the static energy and Floquet spectra [6l |44] . 
This is confirmed numerically in Fig. [ija). It also agrees 
with a broad-band approximation when 17 3> W, using 
the Baker-Campbell-Hausdorff formula, THp = tiHi + 

t2H2 + tit2[H2, Hi] + Y2^1*2[^2-ff2 " ^1-^1, [H2, Hi]] + ■ ■ • . 

The first two terms yield a static Kitaev chain with an 
averaged chemical potential ((/i)) = [ti^i +^2/^2)/^. The 
only terms contributing to the commutator [H2,Hi] are 
AcTj; in Ar^r+i and —ifiay in Ar^r, yielding a Uz term 
in Ar^r+i that contributes to Imw. Physically, Imw in- 
troduces a supercurrent in the chain, which renormalizes 
the spectral gap but, when small, leaves the topological 
phase boundary unchanged ^38j . The same is true for the 
next term shown. We conclude that, when 2 > W, 
there are no Floquet Majorana fermions in the broad- 
band limit. 

It is worth pointing out here that the quasienergy 
gap, Eg, in the Floquet spectrum should not be naively 
used in comparison with other energy scales in the prob- 
lem. For example, even when the average Hamiltonian 
{tiHi +t2H2)/T supports Majorana fermions, eg may be 
vanishingiy small, as illustrated in Fig.jljb) for the static 
limit fii ^ fi2- In this case, the correct energy scale en- 
tering the adiabatic evolution is given by the actual en- 
ergy gap. We shall see eg is relevant for the Floquet sum 
rule, though not for the conductance itself. 

Transport. — Electrons in a driven system do not fol- 
low the usual statistics in a closed system. In the trans- 
port problem we can address this issue by assuming that 
the leads are static and follow the usual Fermi-Dirac 
statistics. The scattering problem between the leads can 
then be formulated by integrating out the leads using 
their Green's function [45l |46] . This procedure adds to 
H„ an imaginary self-energy iT{t) = iEA^^W' where 
T^{t) = 2Im [K^''g^{t)K^] , and is the Green's func- 
tion of lead A. The wire's Green's function is periodic, 
G{t + T,t' + T) = G{t, t'), and satisfies 

[dt - A{t)]G{t, t') - (F * G){t, t') = -i5{t - t'), (3) 

where F * G is the convolution jl^^{t — s)G{s,t')ds. 
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FIG. 2. Differential conductance cr — dl /dV vs. bias V in a 
single-terminal setup for a static system with L — 80, A/w = 
0.6, fi/w — 0.25 and v/w = 7r/50. The (blue) dots are calcu- 
lated from the analytical expression of peak heights. 



Then the steady state (time-averaged) current in lead 
A, = ie {{[H {t) , a^^ a^]J) , can be computed with 
the Green's function. Assuming the leads' density of 
states is constant over the scattering energy range, 
we find an energy-independent 

dJ^/dVx, with bias V\, reads 



= 2ttK^^ p'^K^ and the differential conductance 



27r 



'- n 

+ Tr[e^G(")(..)e^^G(")n^)]/;(-^)} +c^a(^)/;h 

(4) 

where fxiu}) = [I + ef'^"''^^)/''^]"^ is the Fcrmi-Dirac dis- 
tribution, G("n^) ^ tIo I e"*^*e*"''G(i,t-s)dsdt, and 
CaH = E«^A,„Tr[(e+?''^)G(")(c.)e^G(")t(..)] gB]. 
The static case is found by setting G("^(w) — SnfiG{u}). 

For a single lead with a point contact, G vanishes iden- 
tically and the contact matrix is zero except at the con- 
tact site where ^ = v{l — (Jy)/2 with ly — 2ttp\w\'^. Then, 

n 

(5) 

where the Qeh and Ghe are the off-diagonal elements of 
the Nambu-Gorkov Green's function at the contact site, 



I (n+fe) 



(fe)| 



(6) 



with —i5a the self-energy correction to quasienergy e^, 
and and ((^q|, respectively, the right and left Flo- 
quet eigenvectors of the effective (non-Hermitian) Hamil- 
tonian + *r at level a. 

In the weak-contact limit v/w <^ 1, we can employ 
perturbation theory in F. To the leading order, we find 
Wa) = \4>a), (<^a| = {4>a\ (i-S. the Same as eigenvectors 
of H^), and 5a = — ((0a|r|0Q)) [5]. Let us first work 



out the static case. Then, 5a = \v{\u'a\^ + |wqP) with 
Ua and Va evaluated at the contact site. At zero tem- 
perature, limy_j.£;^ '^{y) = o'aL{ ^7^° ) where Ea is an 
energy level of the static system, L{z) = (1 + z^)^^ is the 
Lorentzian, and the peak value. 



2e^ 
27r 



2u'i 



,,c 12 



r,c 12 



(7) 



For the zero-energy Majorana fermion, uq = Vq, so ctq = 
2e^//i in restored units, as is well known [531 HZ]- In 
Fig. [2j we compare this analytical expression with a full 
numerical solution. 

Floquet sum rule. — In the driven system, the peaks at 
V — e{) and V — t.,^ are not quantized even when Floquet 
Majorana fermions are present. This is because energies 
ta + nO. are all connected via the drive force by emission 
and absorption processes. Instead, we find a "Floquet 
sum rule" for the sum of differential conductance at these 



energies 



a{V) = ^cr(t/ + nrj). 



At zero temperature, limy_j.e^ cr(^) = '^aL{^ 
again, a Lorentzian with the peak value. 



(8) 



'-) is, 



2e2 
27r 



WVaf 



(9) 



where ||z|p = X^fc l-^'''^'' particle-hole symmetry 
ll^oll = ||vo|| and IImttII — Wv-kW- Thus, if there is a Floquet 
Majorana fermion at eq and/or e^. 



CTo = cr(eo) 



2e' 



and/or a.^ = cr(e7r) 



2e2 
1^' 



(10) 



respectively. These relations can be generalized for non- 
point contact terms as well. This is our central result. 

The two Floquet Majorana fermions overlap and split 
away from cq or e,r by an amount A that is exponentially 
small in their separation. When \ > v, a also splits 
with the peak values ao and ct^ shifting to = eq i A 
and 1^ = e^r ± A, respectively, each with half the widths 
of the central peak. Therefore, the total weight stays 
the same. This is a general feature: the total weight 
a{V)dV = a{V)dV (x E„ / l^?i,:^(^)Pd.. is con- 
stant at all temperatures. 

We have numerically investigated the Floquet sum-rule 
quantization in the quantum wire. The plots in Fig. [3] 
show the steady differential conductance calculated for 
a two-lead setup with symmetric biases ±V. It is clear 
that, within our numerical precision, ctq and/or 5^ are 
quantized at 2e^/h exactly when eo and/or Floquet 
Majorana fermions appear. Note that the individual 
peaks of cr at y = nfl (for cq Floquet Majorana fermion) 
or V = {2n + l)il/2 (for e-^ Floquet Majorana fermion) 
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FIG. 3. Differential conductances a (top row) and a (bottom row) of the driven system as a function of bias V/r2 in a two- 
terminal setup. The parameters A,^i,/i2 and ti/t2 are as in Fig. [l][a), i//w = 7r/25, and the other parameters are: (a,d) 
L = 40, n/W = 0.37, (b,e) L = 70, Q./W = 0.49, (c,f) L^40,Q/W ^ 0.75. The frequencies are marked by arrows in Fig.jlja). 
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FIG. 4. Differential conductance in a two-terminal setup av- 
eraged over 50 disorder configurations. The parameters are 
as in Fig.|3](a,d) and disorder strength /i^ ~ 0.07VK = 0.19n. 



are not quantized. Indeed, the main contribution is not 
even from n = |46j . The Floquet spectrum is naturally 
reflected in a: The quantized peaks at eq (or e^r) are 
separated from the other peaks by a value of ^ ^ e^^o 
(or Eg.Tr), i.e. the gap in the quasienergy gap separating 
the respective Floquet Majorana fermions from the other 
states. The quasienergy gaps in Fig. [Sjare ~ O.lfi except 
for eg,o - Omn in Fig.[|;e). 

Effects of disorder. — The natural question to answer 
at this point is whether and how a could be measured in 
an actual experiment. Especially important is to be able 
to tell apart a quantized peak from the other features, 
which is complicated if e^^o and e^^^ are small. A pos- 
sible way around is to exploit the topological character 
of the quantization of ctq and cfj^ . Specifically, they must 
be protected against disorder while the other features 



are not. We have studied the effects of disorder numeri- 
cally by adding a static, uncorrelated, random dfir to the 
wire's chemical potential at site r, i.e. (^Mr)dis ~ ^ 
{SfJ.rSnr')d;s ~ t^'d^rr' whcrc is the disorder strength. 
A typical result for the disorder-averaged (cr)^;^ at mod- 
erate disorder is shown in Fig. |4] The central quantized 
peak remains nearly unchanged, while the other peaks 
are suppressed significantly. Note that here /id > Cg^o- 
For stronger disorder the quantized peak is suppressed 
as well pS] . 

Concluding remarks. — In sum, we find that the differ- 
ential conductance summed over periodic drive harmon- 
ics, CT, signals Floquet Majorana fermions with a topo- 
logically protected quantized value 2e^ /h at the Floquet 
Majorana quasienergy. The quantization is robust and 
most prominent in the presence of weak disorder. This 
suggests disorder can be used as a knob to probe Floquet 
Majorana fermions. 

We have studied the finite temperature (r) behavior of 
a and found |46| that for temperature v t < Eg, the 
Majorana peak height falls as ctq.tt ~ 1/t. For t > eg, the 
non-topological peaks start to overlap with the Majorana 
peak and cto.tt saturates. 

In closing, we sketch some problems for future stud- 
ies. Other transport signatures of Floquet Majorana 
fermions, such as noise and heat transport, are interest- 
ing, open problems. Our inclusion of disorder was only as 
a static component of the chemical potential. This is ap- 
propriate if disorder potential is intrinsic to the wire itself 
and is not affected by the external drive. Other disorder 
configurations, e.g. in the contacts or the external drive 



itself, would be interesting to study in future. In particu- 
lar, disorder in the drive spoils the strict periodicity and 
Floquet structure. Quantifying the effects of disorder is 
another interesting problem: the usual way of relating 
the disorder mean-free path to the density of states in 
the static system does not work here, since the periodic 
chemical potential is never in the topological band, so the 
relevant density of states vanishes. Finally, the effects of 
disorder at finite temperature as well as interactions are 
left to future studies. 

We acknowledge useful communications with H. Fer- 
tig and A. Levchenko. This research is supported by 
the College of Arts and Sciences at Indiana University, 
Bloomington. 
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SUPPLEMENTAL MATERIAL 



over the lead states 



Here we sketch the details for the derivation of the 
conductance and Floquet sum rule described in the main 
text. We employ Green's function approach for deriv- 
ing the charge current in the system described by time 
dependent Hamiltonian H^(t) and contact Hamiltonian 
Hr, as discussed in the main text 



Hr, = ^a^''K^j + h.c. (SI) 



The net charge current flowing across the contact A into 
the wire is {h — 1) 



J\t) = ie [H„(t)+H„N\t)] 
^ie {j'' {t)K^^ (t) -h.c.) 



where is the number operator for electrons in lead A. 
One can use the solution of Heisenberg' equation for the 
electron (in the lead) a'*'(i) 



a\t) = r^\t) + / g\t - t')K^j{t')dt\ (S3) 

J ta^ — oo 



where, to is the switching time, g^{t — t') is the Green's 
function matrix of electrons in lead A. In the wide band 
limit the density of states of lead A is constant for 
the relevant energy scales and in the simplest situation 



= {2ir)H^ypl,h{Lo)5{u: - J), (S4) 
= {2^)Hxx'ptJx{^)5{Lo - u'), (S5) 

where /a(w) = 1 - /a(w) = [i + e^'^^'^^^)/^^] is the 
Fermi-Dirac distribution of the lead A with bias Vx and 
temperature r^. (The Boltzmann constant = !■) 
For the Majorana operator 7(t), the integration of the 
Heisenberg equation can be complicated by the time de- 
pendence of A{t). In a periodically driven system it is 
obtained in terms of the Floquet Green's function, as we 
discuss later. 

Static System 

If the time dependence of A{t) is trivial, one can inte- 
grate the Heisenberg's equation for the Majorana opera- 



(S2) tor with using Eq. (S3 1. In the Fourier space. 



7(c.) = G{u:)h{uj), 

where G(w) is the Green's function defined from 

G-'^{uj) =uj -iA-iV{uj), 
with the self energy 

iT{uj) ^ \K^^g^{u)K^ - K^^ g^*{-uj)K^ 

and 



(S6) 

(S7) 
(S8) 

(S9) 



g (u)) = —inp . The noise term ri {t) — ig {t — to)a {to) Using Eq. (S6| and Eq. (S3) in the current equation 



obeys the fluctuation-dissipation relation after averaging Eq. ( S2 ) and averaging over the lead states one obtains [T] 



= / dLo{Tv[eG{u;)e'G\oj)\ [fx>{oj) - fx{uj)] + [eGioj)^" GHlo)] [1 - fx'{-uj) - fx{u)]} , (SIO) 



where S,^ ~ 2ttK'^^p^K^ and = J {J^{t))&t. The The second term vanishes if the system has a single lead, 
conductance is obtained by taking the derivative with 
respect to the bias Vx, a 



e 
'2tt 



Auj Tr 



eG{uj)e' G^Lo) 



x[/;(-c.)+/;(^)] 



(S12) 



e 

2^ 



dwTr 



eG{u:)e'G\u:) [/;(-^)+/^(^)] 



27r 



X'^X 



J2 /d^Tr eG{Lu){e' +e")G\Lo) fiiu) 



Periodically driven system 



In a driven system, the formulation of current by non- 

(SIX) 

equilibrium Floquet Green's function has been discussed 
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before [2H3] and we follow the same strategy. Integrating 
the Heisenberg equation for the Majorana operator in a 
driven system gives |5] 

^--^Ait) 



-/{t)~il d.sT{s)-f{t-s)^h{t), 



(S13) 



with self energy iV (s) = 2iJ2x^^[K^''g^{s)K^] and 
hit) = 2iJ2x'^^[K^'^V^it)]. The Green's function of 
this inhomogeneous equation G{t,t') satisfies 



G{t,t')~i I dsT{s)G{t~ s,t')ds 



= 6{t-t'). 



(S14) 



For a periodic drive with period T — 27r/J7, the Floquet 
Green's function is also periodic over the same period 



G(t + T,t' + T) = G{t,t'). One can introduce the Fourier 
transform, 



G{t,t') - E / Gik)(^^^^~^u.it~t')^-^knt^ (g^5) 



The Majorana operator is solved in terms of the Green's 
function 



1^ ^ 



— iujt^—ikilt 



G^^\u^)h(uS). (S16) 



Using these expressions in the current formula Eq. (S2|, 
the the steady state current = ((J^)) = 

/o {J^i^)) 'i-t/T in the wide band limit is found to be 



e 
Tr 



5: /dc. ^ (Tr [G(')t(..) {e'+e") G(H-^)e] hH 

"^^G(')Me''GWt(a;)l /^,(^) _ Tr k^G«(c.)C^"G(')t(^)l [i _ fy{-u:)] I (S17) 



which gives the conductance 
'2tt 



a 



A 



j duj ^Tr[e^GWHe''G(W(, 



e' 
2tt 



J2 fdu; ^Tr [e^G«(^)(^^' +^^'')G(')t(,,)] /;(^). (S18) 



Finally, similar to the static problem, for a single lead, 

o2 



"2^ 



^ / d^Tr[e^G«(^)e^'GWt(, 
I 

x[/;(-^) + /;m]. 



r 



with vanishing positive e. 

For a single lead with point contact, the matrix ^ 



v{\ — ay) 12. By unitary transforming the Eq. (S12) with 
U we get 



(S19) 



27r 



dc. |f;ehHl'[/U-'^) + /;M], (S21) 



Relation to Nambu-Gorkov Green's function 

The Majorana representation is related to the Nambu 
spinor (c^,cj^) = "ilil by the unitary transformation 7,. = 

mbr with u — A= ( ^ . ^. ] 

V2 \ —I I I 



where tJch is the electron-hole component of the Green's 
function evaluated at the contact site 



(S22) 



for each site r. And the 



For the periodically driven system the relevant Green's 



Bogoliubov-de Gennes (BdG) Hamiltonian is related as function is the Floquet Green's function in Nambu spinor 
ffsdG = iW'AU, where U = u (B u (B ■ ■ ■ (S u (L times). 
The Nambu-Gorkov Green's function in a static system 
is defined as 



BdG 



^e)- 



(S20) 



basis, which is defined as 5 



(fc)i 



(S23) 
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where \ipa^) and {'fi''a ^\ are the discrete Fourier compo- 
nent of the time periodic function 'Floquet states' \(pa{t)) 
and {(fa (t) I ) which are right and left eigenstates of the ef- 
fective Floquet Hamiltonian i/off = H„ + iF — i ^ with 
eigenvalues [Floquet energies) — iSa- 

HeS l^a(t)) = (fa " iS^) |</?a(t)) , (S24) 
((^a(i)l HeS = (ea - iSa) {'Pait) \ ■ 

\ipa(t)) lives in an extended (with the time axis) Hilbert 
space [6] , where the expectation of any operator O in the 
state is defined with the time average 



k 



(S25) 

In terms of the Green's functions Q^^\ one can express 
the conductance for a periodically driven system with 



single lead and contact matrix ^ 
the static case, 



/^(l — (Tj,)/2 similar to 



(S26) 
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FIG. SI. Contributions to Floquet sum rule from different 
n in a(V) = En + '^)- Here a'") = a{nQ. + V). The 

parameters are as in Fig. 3 (c,f) of the main text. 



Now, the eigenstate of the time dependent Hamilto- 
nian i/w, |'0q(O) ^re related to the Floquet state \(pa{t)) 
by the relation \ipa{t)) = e-*('°~"'°)* For 

and tt/T Floquet energy, |V'o(i)) = I'Poit)) and 

\^p^{t)) = e~"*/^e-''-* |(/?,r(i))- Considering the contact 
term perturbatively, in the leading order one can neglect 
the modification to the wavefunction to have \ipo{t)) « 
\ipo{t)) and « e"*''*/'^ Iv-Kit)). Using the particle- 

hole symmetry at Floquet energies 0, tt/T, one can 
express the eigenstates in the Bogoliubov quasiparticle 
representation IV'ol*)) = (uoit) , ul{t)y and |V-'7r(0) = 
{uTr{t),u%{t)y . This allows us to write \<fi{t)) in the 
Fourier space 





(S27) 



The leading perturbation to the self-energy part can be 
computed in the extended Hilbert space [5] with 



F = -(e + e)/2 = -, 



1 
1 



@ p® . ..® q , (S28) 

L — 1 times 



Sa^- {{(pa{t)\r = -(iKip + \K\\'), (S29) 

where u^'^^ and Va ''^ are evaluated at the contact site 
and ||z|p = Xlfc k^'^^P- We also note that 



c(fc-l-i) c(fe) 



IKIP + IKIP 



y — -I- 
V-ea 



5a 



(S30) 



where ^(z) = (1 -I- z^) ^ is the Lorcntzian. The peak 
value at zero temperature is 



cTq = ^lim ^ a{V + fcJl) = 



2e^ 
2^ 



2IKII IKI 



|K|P + ||«c||2 



(S31) 



This is our Eq. (9) in the main text. 



Numerical Results 



In Fig. SI we plot u(V) = ^„ aiy + nil) for system 
with two symmetrically biased leads {Vi = — V2 — —V) 
hosting a eg Floquet Majorana fermion, which shows that 
the peaks of individual components a{V + nfl) are not 
quantized, but the sum a shows a quantized peak. 

In Fig. |S2[ we compare the differential conductance 
a for smaller and larger than optimal disorder strength. 
For small disorder, the non-topological peaks do not get 
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FIG. S2. Differential conductance in a two-terminal setup averaged over 50 disorder configurations, where disorder strengths 
(a) fid = 0.05W^ = 0.13O and (b) fid = O.llVK = O.SOil which are smaller and larger than the disorder strength used in the 
main text. Other parameters are as in Fig. |S1[ except here u/w — tv/17. We show the contribution to the sum, a-{V) = E, 
from different a''"\V) = aiy + nfJ). 



suppressed enough compared to the topological peak, 
whereas for large disorder the topological peak can get 
suppressed. 

Finally, in Fig. [S3] we plot the temperature dependence 
of (t(0) with temperature (r). As r becomes comparable 
to the Floquet gap e^,, the peak becomes flat. 
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FIG. S3. Temperature (r) dependence of the peak height of 
conductance 5(0) = ctq for Floquet Majorana. The param- 
eters used are same as that of Fig. |S1| The Floquet gap is 
e„ = 0.15!r2. 



